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A  iine-dimensmnal  oceanic  planetary  boundary  layer  model  is  used  to  Investigate  the  response  of  the 
upper  ocean  to  the  atmospheric  conditions  which  are  predicted  to  develop  following  a  hypothetical 
nuclear  exchange  The  ocean  model  is  driven  by  the  surface  heat  and  momentum  fluxes  predicted  by  an 
atmospheric  general  circulation  model  following  a  summertime  injection  of  1.5  x  10’ *  g  of  smoke  from 
postwar  tires  over  Europe,  Asia,  and  North  America.  Although  the  specific  response  of  the  upper  ocean  is 
highly  dependent  on  the  geographic  location,  the  mid-latitude  summertime  mixed  layer  typically  cools  .1 
to  5  C  and  deepens  25  m  during  the  first  .10  days  following  the  smoke  injection.  Moreover,  a  large 
fraction  of  this  response  is  found  to  take  place  during  a  short  2-  to  .1-day  period  of  very  intense  winds 
and  falling  air  temperatures,  which  occurs  during  the  first  week  or  two  after  the  smoke  injection. 


Introdix  TION 

In  the  past  5  years  much  attention  has  been  given  to  the 
potential  climatic  consequences  of  a  major  nuclear  exchange 
between  the  superpowers.  One  result  of  this  attention  has  been 
the  ilevelopment  of  the  “nuclear  winter"  hypothesis  that 
tremendously  large  amounts  of  smoke  injected  into  the  atmo¬ 
sphere  from  the  tires  from  hundreds  of  nuclear  detonations 
during  a  major  nuclear  exchange  would  block  enough  sunlight 
that  the  interiors  of  the  northern  hemisphere's  continents 
would  cool  '  10  .10  C  for  periods  of  the  order  of  a  month. 
This  hypothesis  has  developed  as  a  result  of  the  initial  studies 
by  Crtiizen  anil  Birks  [19X2],  who  calculated  the  amount  of 
smoke  that  is  produced  by  a  large  number  of  massive  fires, 
and  by  lurid  pi  ai  [19X3],  who  investigated  the  one- 
dimcnsional  radiative-convective  response  of  the  atmosphere 
to  certain  smoke  aerosol  injection  scenarios.  The  major  con¬ 
clusions  of  Turco  et  al.  have  been  substantiated  and  further 
refined  by  two-dimensional  IMacCriuken.  1983]  and  three- 
dimensional  \ Ali’ksiuulroi  anj  Sienihikov.  1983;  C'oix  v  el  al, 
1984;  Cess  el  al.  1985]  atmospheric  simulations  with  fixed 
distributions  of  smoke  and  by  two-dimensional  [^Haherle  et  al. 
1985  I  and  three-dimensional  [ MaeCrai  ken  anil  ii'alton.  1984; 
Stem  hikin .  1985;  Malone  el  al.  1986;  Thompson.  1985]  atmo¬ 
spheric  simulations  in  which  smoke  is  transported  and  re¬ 
moved  by  the  evolving  atmospheric  circulation.  In  the  three- 
dimensional  simulations,  both  the  magnitude  and  the  geo¬ 
graphical  distribution  of  surface  air  temperature  reductions 
are  modified  by  the  vast  (assumed  infinite  in  most  simulations! 
heat  capacity  of  the  underlying  oceans. 

Since  the  oceans  have  been  found  to  play  such  an  important 
role  in  the  atmospheric  simulations  of  nuclear  winter,  it  is 
surprising  that  with  the  exception  of  the  Soviet  work,  no 
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serious  attempt  has  yet  been  made  to  investigate  the  ocean's 
respon.se  to  such  nuclear  war  scenarios.  Rohoek  [1984]  exam¬ 
ined  the  nuclear  winter  hypothesis  using  a  seasonal  energy 
balance  climate  model  based  on  that  of  Sellers  [1973]  coupled 
to  an  ocean  model  with  a  constant  (in  time)  mixed  layer  depth. 
In  response  to  plausible  nuclear  winter  scenarios,  the  coupled 
model  predicted  sea  surface  temperature  (SST)  decreases 
'll  C  in  100  days  in  mid-latitudes.  In  the  present  study  we 
use  a  physically  realistic  one-dimensional  ocean  planetary 
boundary  layer  model  to  investigate  the  short-term  conse¬ 
quences  in  the  ocean  of  a  three-dimensional  atmospheric  gen¬ 
eral  circulation  model's  response  to  smoke  generated  by  a 
full-scale  nuclear  exchange.  We  consider  only  a  northern 
hemisphere  summer  case  (July)  because  previous  atmospheric 
model  simulations  [rutx'v  ei  al.  1984]  indicate  that  the  atmo¬ 
sphere's  response  to  a  nuclear  exchange  is  much  stronger  in 
summer,  when  more  solar  radiation  is  available  for  absorption 
by  the  smoke,  than  in  winter. 

We  have  focused  our  study  on  the  short-term  respon.se  of 
the  ocean  becau.se  of  the  large  uncertainties  associated  with 
estimating  the  potential  response  of  the  atmosphere  to  a  full- 
scale  nuclear  exchange.  Thc.sc  uncertainties  fall  into  two  major 
categories  [Berber.  1986].  The  first  concerns  the  amount,  dis¬ 
tribution.  and  physical  properties  (stye,  shape,  and  blackness) 
of  the  smoke  produced  during  and  immediately  after  a  nuclear 
exchange.  These  uncertainties  are  due  to  our  limited  knowl¬ 
edge  of  the  specific  details  of  the  nuclear  exchange,  the  avail¬ 
ability  of  combustible  material,  the  thermodynamics  and 
chemistry  of  large  fires,  and  the  properties  of  particulate 
matter  in  the  smoke  plumes.  The  second  category  concerns 
the  ensuing  elfect  of  the  smoke  layer  on  the  atmosphere  and 
the  elfect  of  the  atmosphere  in  dispersing  or  removing  the 
smoke  as  well  as  changing  the  optical  properties  of  smoke 
particles  by  aggregation  or  chemical  reactions.  As  a  result, 
reliable  quantitative  estimates  of  the  atmospheric  response  to 
a  nuclear  exchange  are  probably  limited  to  less  than  I  month. 
For  such  short  time  scales  and  for  the  rather  strong  atmo- 
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I  ig  I  The  milial  dislrihulion  of  smoke  extinction  optical  depth 
Used  in  the  atmospheric  nuclear  winter  simulation  [from  Glum  ei  ui, 
lds5J,  Dots  show  the  locations  of  the  ocean  numerical  experiments, 
and  the  number  beside  each  dot  gives  the  dilTerences  in  SST  (in  de¬ 
grees  C'elsiusi  after  ,^0  days  between  the  control  and  nuclear  winter 
experiments.  At  all  locations  the  SST  was  lower  after  the  nuclear 
winter  experiment. 

spheric  forcing  that  is  expected,  the  response  of  the  upper 
ocean  clearly  will  be  dominated  by  the  local,  i.e.,  the  one- 
dimensional,  forced  re,spon.se. 

Tuf  Moofx  and  thk  Nlmfrical  EXPKRI.MhSrS 

In  this  study  we  use  the  one-dimensional  ocean  planetary 
boundary  layer  model  of  Garwood  [1977],  The  model  uses  the 
turbulent  kinetic  energy  budget  to  describe  the  production, 
dissipation,  and  buoyant  damping  of  turbulent  kinetic  energy- 
in  the  upper  ocean.  Separate  equations  for  the  vertical  and 
horizontal  components  of  turbulent  kinetic  energy  are  used, 
and  closure  is  achieved  by  parameterizing  the  second-  and 
third-order  terms  in  the  turbulent  kinetic  energy  equations. 
While  the  Garwood  model  has  a  number  of  unique  features 
which  make  it  somewhat  different  from  other  bulk  or  profile 
models  of  the  upper  ocean,  other  modern  mixed  layer  models 
are  expected  to  produce  qualitatively  similar  results  [Marlin. 
1985]. 

Two  types  of  numerical  experiment.s.  "control”  and  "nuclear 
winter."  respectively,  were  carried  out  at  six  different  locations 
in  the  North  Atlantic  and  North  Pacific  oceans.  In  the  control 
experiments,  the  ocean  model  was  driven  by  the  local  wind 
stress  and  heat  fluxes  at  the  sea  surface  taken  from  a  "normal" 
or  "typical"  atmospheric  general  circulation  model  .simulation 
of  the  "perpetual"  month  of  July.  In  the  nuclear  winter  experi¬ 
ments.  the  surface  heat  fluxes  and  surface  stresses  were  taken 
from  an  otherwise  identical  atmospheric  model  simulation 
except  that  a  total  of  1.5  x  10''*  g  of  smoke  was  injected  into 
the  air  tiver  Europe.  Asia,  and  North  America,  The.se  atmo¬ 
spheric  simulations,  which  made  use  of  the  Oregon  State  Uni¬ 
versity  two-level  tropospheric  general  circulation  model  [Sch- 
lesincjcr  and  Gales.  1980;  Ohan  el  al.  1982],  have  been  de¬ 
scribed  in  some  detail  by  S.  J.  Ghan  ct  al.  (The  climatic  re¬ 
sponse  to  large  atmospheric  smoke  injections:  .Sensitivity 
studies  with  a  tropospheric  general  circulation  model,  submit¬ 
ted  to  Journal  oj  (leophysieal  Re.search.  1986).  In  both  atmo¬ 
spheric  simulations,  control  and  nuclear  winter,  the  sea  surface 
temperatures  were  prescribed  to  be  July  climatological  values. 
Tor  these  studiev  the  atmospheric  model  was  mndib'’d  by  the 
incorporation  of  a  deita-Eddington  formulation  for  solar  radi¬ 
ation  ff’ess  el  al..  1985],  and  in  the  nuclear  winter  simulation 
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it  was  coupled  to  a  Lagrangian  trace  species  transport  model 
known  as  GRANTOUR  (J.  J.  Walton  et  al.,  A  global-scale 
Lagrangian  trace  species  transport  model,  sunmitled  to  Jour¬ 
nal  of  Geophysical  Research.  1986).  By  means  of  GRAN- 
TOLi'R,  smoke  particles,  which  were  initially  distributed  uni¬ 
formly  with  height,  were  subjected  to  dry  deposition  and  co¬ 
agulation  as  well  as  precipitation  scavenging  and  large-scale 
advective  transport.  The  particular  nuclear  winter  experiment 
considered  here  assumes  an  injection  of  1.5  x  lO'-*  g  of  smoke 
with  a  refractive  index  of  in  =  1,75  —  O..Ti  and  a  lognormal  size 
distribution  with  a  number  mode  radius  of  0. 1  pm  and  a 
standard  deviation  of  2.0.  Figure  1  shows  the  initial  distri¬ 
bution  of  the  smoke  as  represented  by  the  extinction  optical 
depth  [Turco  ei  al..  1983],  Values  greater  than  50  are  found 
over  both  continents.  By  the  end  of  the  30-day  simulation,  the 
smoke  concentration  becomes  relatively  uniform  over  the 
northern  hemisphere,  w  ith  values  of  optical  depth  in  the  range 
from  0.5  to  2.0  (S.  J.  Ghan  et  al..  The  climatic  response  to 
large  atmospheric  smoke  injections:  Sensitivity  studies  with  a 
tropospheric  general  circulation  model,  submitted  to  Journal 
of  Geophysical  Research.  1986).  Figure  I  also  shows  the  six 
locations  that  were  used  in  the  ocean  experiments.  They  were 
chosen  to  represent  mid-latitude,  subtropical,  and  tropical  re¬ 
gions  and  to  represent  a  variety  of  distances  from  the  initial 
sources  of  smoke. 

Since  no  present-day  atmospheric  general  circulation  model 
simuia'cs  the  real  atmosphere  perfectly,  it  was  not  surprising 
to  find  that  the  mean  wind  stress  and  surface  heat  fluxes  com¬ 
puted  in  the  atmospheric  control  simulation  differed  some¬ 
what  from  climatology.  Therefore,  to  provide  more  realistic 
atmxvspheric  forcing  to  the  ocean  in  the  control  experiment 
and  presumably  also  in  the  nuclear  winter  experiment,  the 
wind  stress  and  heat  fluxes  from  the  atmospheric  simulations 
were  calibrated  to  climiatology  as  follows.  The  mean  surface 
wind  stress  r,  downward  solar  radiation  Q,,  long-wave  radi¬ 
ation  Q^.  and  latent  ((>,.)  and  sensible  (Q^)  heat  fluxes  from  the 
atmospheric  model  control  simulation  were  compared  with 
the  corresponding  July  climatological  values  of  Eshensen  and 
Kushnir  [1981],  If  f  denotes  one  of  these  quantities  (the  wind 
stress  or  a  heat  flux  component),  a  calibration  constant. 

Gf  =  F„  F,  (1) 

was  computed  for  each  F  at  each  of  the  ocean  locations  shown 
in  Figure  1.  In  (I|.  the  overbar  represents  an  average  over  the 
month  of  July,  and  the  subscripts  "o"  and  “v"  denote  observed 
[ Eshensen  and  Kushnir.  1981]  and  simulated  (control  experi¬ 
ment).  respectively.  The  resulting  values  of  Cf  are  shown  in 
T  able  1.  Many  of  the  values  are  near  unity,  indicating  that  the 
mean  surface  fluxes  from  the  July  control  simulation  agree 
quite  well  with  July  climatology.  The  anomalous  value  of 
r,  =  -5.  associated  with  F  =  (J.,  at  point  3,  is  caused  by  the 
occurrence  of  very  small  values  of  Q^.  at  that  location.  A  cali- 
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lirateJ  Iiiiil'  of  ihe  Miilaee  tlu\cs  at  cacli  ItK'ation  was 

eon^l^Ue'lcd  h\  niultipl\ ing  the  (i-tuniiK  \ allies  of  r.  O^.  Q^.  O  . 
and  (Ji..  produced  h\  the  atitiospherie  iiiotlel  siimilations.  hy 
the  appropnale  ealihratioii  constant  in  Talile  1.  Note  that  as  a 
consei.|uence,  hoih  the  mean  and  the  standard  deviation  of  the 
Muses  arc  moditieil  h>  the  same  factor,  t  , .  This  procedure 
eu.iiaiuces  (iiat  (h  ■  time  scries  of  each  eontponeiit  of  the  siir- 
t.iec  forcing  that  drives  the  ocean  model  in  the  control  experi¬ 
ment  has  as  its  mean  the  corresponding  observed  cliiiiatologi- 
c.d  me. in  for  .Iiilv.  This  constraint  on  the  mean  vvf  the  atmo¬ 
spheric  fvMcing  terms  is  needed  to  assure  a  realistic  control 
st.ite  in  the  ocean  H>  using  the  same  calibration  factors  in 
both  the  control  and  the  nuclear  winter  expermient.  we  cssen- 
liallv  assume  that  the  biases  in  the  atmospheric  model  lUixcs 
are  unaDeeted  by  nuclear  winter. 

Table  ’  shows  a  summary  of  the  calibrated  surface  Iluxes 
that  were  Used  in  the  ocean  experiments.  It  can  be  seen  that 
the  .iiniosphene  nuclear  winter  experiment  produced  signill- 
c.inl  changes  in  the  surface  wind  s'ress  and  heal  Iluxes  at  most 
ol  the  geographic.il  locations  studied.  More  speeilieally.  the 
.iver.ige  wind  stress  increased  lalmost  doubled!  at  four  lo- 
c, molls,  while  the  net  heal  Ilux  decreased  at  all  six  locations. 
I  he  l.irgesi  deere.tse  in  net  heat  Ilux  occurred  at  the  three 
oce.in  loc.itioiis  ipoinis  2,  T  and  bl  lying  directly  downwind  of 
ihe  sinohc  sources,  {)n  the  other  hand,  much  smaller  decreases 
occuricil  m  the  e.istern  Nortli  Paeilie  Ocean  (points  1.  4.  aiul 
■')  I  hese  dcsre.ises  in  the  net  surface  heat  Ilux  were  largely 
due  to  ilecre.ised  sol.ir  r.idialion  caused  by  the  direct  absorp¬ 
tion  ol  sol.ir  r.Kli.ilion  by  the  smoke.  However,  at  point  6. 
which  Is  loc.ited  m  the  hast  (  hiria  Sea.  the  decreased  net  heat 
Ilux  Is  pi  im. inly  ilue  to  .i  very  hirge.  ami  to  some  extent  unre- 
.ilislic  isee  below  i.  increase  in  the  surface  evaporation. 

\i  c.ich  loe.ilion.  the  oce.in  model  was  driven  by  the  cali- 
isr.iled  V, lines  o|  ihe  surf, ice  Iluxes  computed  every  b  hours 
tioin  the  .ilmosphcrie  simulations.  The  precipitation  rate  from 
the  sinuil.il ions  w.is  .iKo  used  to  compute  the  surface  salinity 
lluv  needed  by  Ihe  oce.in  moilel.  However,  for  brevity  and 
bec.iuse  the  response  of  the  salinity  in  the  nuclear  winter  ex- 
periments  w.is  generally  small,  neither  precipitation  nor  salini¬ 
ty  will  be  discussed  m  this  note. 

I  he  initial  conilitions  for  the  ocean  experiments  consisted  of 
temperature  ami  salinity  proliles  from  the  surface  to  a  viepth 
of  200  m  Ihe  same  initial  conditions  were  usevi  in  both  the 
control  experiment  .ind  the  nuclear  winter  experiment  Ihe 
initial  conditions  were  based  on  climatological  protiles  of  tem¬ 
per, iture  ami  s.ilimty  for  iuly  at  each  location,  obtained  from 
fleet  Numerical  Oce.inography  Center  (I  NOC)  in  Monterey. 
(  .ilifornial  Model  sensitivity  tests  indicate  that  the  response 


of  the  upper  oeean  in  these  experiments  is  not  sensitive  to  the 
particular  initial  conditions  used  as  long  as  they  are  reason¬ 
ably  representative  of  Iuly  conditions  at  the  location  in 
question. 

Oc  t  \\  Rise  I  It; 

The  response  I'f  the  sea  surface  temperature  (.S.STI  and 
mixed  layer  depth  iMI.Dl  in  the  control  and  nuclear  winter 
experiments  is  shown  in  Figures  2  and  .V  Looking  lirst  at  the 
SSL  (Figure  2|.  we  see  that  at  each  location  studied  the  tem¬ 
perature  is  lower  at  the  end  of  the  nuclear  winter  experiment 
than  It  is  at  the  end  of  the  control  experiment.  The  drop  in 
SSI  due  to  the  nuclear  exchange  is  typically  .f  .“v  C.  although 
it  IS  less  than  I  C  at  two  locations  and  more  than  S  C  at 
another.  Where  a  signiticani  dillerenee  develops  between  the 
control  and  nuclear  winter  SS  I's,  the  dillerenee  develops  quite 
suddenly,  usually  over  a  period  of  only  a  few  days.  At  two 
such  locations  (,s0  N.  l.xO  W  and  .xO  N.  150  FI.  the  dilTerence 
devehvps  because  of  a  sudden  anomalous  cooling  of  the  upper 
ocean  which  takes  place  during  the  nuclear  winter  experiment, 
while  at  two  other  locations  the  dilTerence  develops  because  a 
short  periovi  of  strt'ng  warming  which  occurs  in  the  control 
experiment  does  mvt  occur  in  the  nuclear  winter  experiment. 
This  dilTerence  between  the  control  and  nuclear  winter  SSTs 
begins  to  develop  at  a  dilTercnt  time  at  each  location,  depend¬ 
ing  on  the  distance  upw  ind  to  the  nearest  source  of  smoke.  In 
mid-latitudes,  where  many  of  the  experimental  points  are  lo¬ 
cated.  the  smoke  is  transported  generally  eastward.  Thus  the 
smoke-induced  cooling  begins  after  only  T  days  in  the  western 
North  Atlantic  olT  the  coast  of  Newfoundland  (5t)  N.  50  W). 
but  it  docs  not  start  until  about  day  12-15  in  the  eastern 
North  Pacilic  (50  N.  150  W),  The  locations  farther  south  in 
the  North  Pacilic  t.tO  N.  150  W  and  10  N.  150  W)  arc  not 
only  f.ir  viownvvind  but  also  far  to  the  south  of  the  snuike 
source.  As  a  result,  the  response  of  the  upper  ocean  to  the 
smoke  injection  is  very  small  at  these  locations. 

At  most  of  the  locations  shown  in  Figure  2.  the  dilTcrences 
between  the  control  and  nuclear  winter  SS  T  s  bece  me  rela¬ 
tively  constant  in  lime  after  about  2  or  .2  weeks.  An  important 
exception  to  this  behavior  occurs  at  the  pviint  located  in  (he 
l  ast  C  hina  Sea  t.fO  N.  125  Fl.  where  the  SST  in  the  nuclear 
winter  experiment  is  still  falling  in  relation  to  the  control  ex¬ 
periment  after  TO  days.  As  was  noted  above  in  connection  with 
the  surface  Iluxes  (Table  2l,  the  nuclear  exchange  has  a  rather 
unique  elTecl  on  the  surface  heat  Iluxes  at  this  location.  Here, 
in  the  atmospheric  luielear  winter  simulation,  the  atmosphere 
extracts  extremely  large  amounts  of  sensible  and  latent  heat 
fmm  (he  undeilying  oeean.  V\e  t.in  see  from  Table  2  that  the 
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an  important  first  step  toward  the  ultimate  goal  of  determin¬ 
ing  the  true  response  of  the  coupled  ocean-atmosphere  system. 

Turning  now  to  the  MLD  (Figure  3),  we  see  that  at  a  ma¬ 
jority  of  the  locations  studied,  the  MLD  generally  becomes 
deeper  during  the  nuclear  winter  experiment  than  during  the 
control  experiment.  As  shown  here,  the  MLD  refers  to  the 
turbulent  boundary  layer  depth,  a  prognostic  variable  in  the 
ocean  model.  Being  the  instantaneous  depth  of  penetration  of 
turbulent  mixing  in  the  upper  ocean,  the  MLD  will  not  neces¬ 
sarily  correspond  to  the  conventional  definition  of  mixed  layer 
depth,  namely,  the  depth  at  which  the  temperature  drops  to 
the  value  0.2  C  below  the  SST,  although  there  is  a  strong 
similarity  between  the  two.  The  distinction  is  especially  rele¬ 
vant  during  those  periods,  usually  of  short  duration,  when  the 
MLD  has  shallowed.  In  these  cases,  the  temperature  drop 
below  the  MLD  can  be  very  small.  Nevertheless,  at  the  three 
high-latitude  points  studied  (i.e.,  at  50  N),  a  noticible  deep¬ 
ening  of  the  oceanic  boundary  layer  occurs  in  the  nuclear 
winter  experiment  at  the  time  when  the  corresponding  SST 
(Figure  2)  first  begins  to  differ  from  the  control  SST.  Off  New¬ 
foundland  (50  N,  50  W),  however,  the  MLD  slowly  recovers 
from  its  sudden  response  to  nuclear  winter,  and  by  the  end  of 
the  experiment  it  is  very  similar  to  that  in  the  control  experi¬ 
ment.  At  the  East  China  Sea  location,  the  MLD  in  the  nuclear 
winter  experiment  does  not  shallow  on  day  10  as  it  does  in  the 
control  experiment,  and  this  is  when  the  two  SSTs  (control 
and  nuclear  winter  (Figure  2))  suddenly  begin  to  depart  from 
each  other. 


Time  (days) 

Fig  2.  Time  evolution  of  the  SST  at  each  ocean  location  in  the 
control  (solid  linel  and  nuclear  winter  (dashed)  experiment. 


increased  sensible  plus  latent  heat  flux  at  this  location  ac¬ 
counts  for  259  W  m“^  of  the  342  W  m"^  decrease  in  the  net 
sutface  heat  flux  which  takes  place  during  the  nuclear  winter 
experiment.  Such  extreme  surface  fluxes,  which  do  not  occur 
at  any  of  the  other  locations  studied,  are  caused  in  part  by 
strong  winds  (see  r  values  in  Table  2)  and  very  cold,  dry  air 
which  develops  upstream  over  the  nearby  Asian  continent 
during  the  atmospheric  nuclear  winter  simulation.  Another 
reason  for  the  large  heat  loss  at  the  surface  is  the  fact  that  the 
SSTs  are  quite  high  in  the  East  China  Sea  in  July,  and  these 
high  -SSTs  are  held  constant  throughout  the  atmospheric  nu¬ 
clear  winter  simulation.  Because  of  this,  the  sensible  and  latent 
heat  fluxes  at  this  location  and  our  ocean  model’s  predicted 
large  response  may  be  somewhat  unrealistic.  In  a  syn¬ 
chronously  coupled  atmosphere-ocean  model  simulation  of 
nuclear  winter,  in  which  the  SSTs  are  allowed  to  change  in 
response  to  the  computed  surface  fluxes,  these  fluxes  would 
certainly  tend  to  decrease  as  the  SST  decreases.  Such  de¬ 
creased  fluxes  would  act  to  reduce  the  SST  changes,  but  they 
could  also  alter  the  subsequent  atmospheric  circulation.  It  is 
therefore  not  possible,  using  our  present  results,  to  estimate 
with  any  certainty  what  the  ocean’s  response  in  such  a  cou¬ 
pled  nuclear  winter  experiment  would  be.  In  the  absence  of 
additional  information  or  further  model  tests,  we  feel  it  is  safe 
to  assume  that  the  predicted  response  of  the  SST  at  30’N, 
125  F.  in  Figure  2  represents  an  upper  bound  for  the  change 
that  might  actually  occur  at  that  location  during  a  nuclear 
winter.  Knowing  the  upper  bound  on  the  ocean’s  response, 
and  knowing  some  of  the  ocean  processes  that  are  involved,  is 
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Fig.  .1.  Same  as  Figure  2  except  for  mixed  layer  depth. 
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hig  4  Time  exoliition  of  ihe  calibraied  lal  net  surface  heat  flux  and  (hi  surface  wind  stress  r  at  50  N.  150  E  used  in 
the  control  experiment  (solid  lines)  and  Ihe  nuclear  winter  experiment  (dashed). 


By  examining  Table  2  and  f  igures  2  and  ,7,  it  is  possible  to 
determine  for  each  geographical  location  whether  the  model- 
simulated  SST  decrease  during  the  nuclear  winter  experiment 
was  due  to  enhanced  air-sea  cooling  or  entrainment.  For  ex¬ 
ample.  at  the  first  location  (50  N.  150  W).  the  net  heat  flux 
into  the  ocean  was  55  W  m“*  less  during  the  nuclear  winter 
experiment  than  during  the  control  experiment  (Table  2).  If 
the  average  MLD  at  this  location  during  nuclear  winter  is 
taken  to  be  50  m  (Figure  .7),  this  decrease  in  surface  heat  flux 
cannot  account  for  more  than  a  1  C  difference  in  SST  at  the 
end  of  the  month.  This  implies  that  well  over  half  of  the  2.8  C 
decrease  in  SST  in  the  nuclear  winter  experiment  (Figure  2) 
was  caused  by  entrainment.  This  interpretation  is  also  consis¬ 
tent  with  the  much  greater  wind  stress,  an  important  factor  in 
producing  entrainment,  experienced  during  nuclear  winter  at 
this  location  (Table  2).  In  the  same  way.  we  find  that  more 
than  half  of  the  decrease  in  SST  at  the  second  geographical 
location  (50  N.  I50  El  is  due  to  a  decrease  in  the  net  surface 
heat  flux,  with  the  rest  due  to  enhanced  entrainment  (see 
higure  5/;.  and  additional  discussion  below).  At  the  third  lo¬ 
cation  (51)  N.  50  VV).  where  there  was  little  change  in  the  wind 
stress  and  MLD  due  to  nuclear  winter,  all  of  the  SST  decrease 
is  attributed  to  the  rather  large  decrease  (126  W  m  ^)  in  the 
net  surface  heat  flux.  The  only  other  location  where  the  SST 
response  was  significant  is  the  F.ast  China  Sea  point  (70  N. 
125  E).  and  there  most  of  the  SST  decrease  is  accounted  for 
simply  by  the  very  large  change  in  the  net  surface  heat  flux 
caused  by  the  extraordinary  sensible  and  latent  heat  fluxes 
noted  above. 

We  now  examine  the  response  of  the  upper  ocean  to  nuclear 
winter  in  greater  detail  hy  describing  the  time  evolution  of  the 
vertical  thermal  structure  and  its  relation  to  the  atmospheric 
forcing  at  one  of  the  experimental  ocean  locations.  The  point 
at  50  N.  150  Ei  is  chosen  for  this  purpose  because  the  atmo¬ 
spheric  forcing  during  the  control  and  nuclear  winter  experi¬ 
ments  at  this  location  is  fairly  representative  of  the  other  lo¬ 
cations  as  well  This  point,  located  approximately  500  km  west 
of  the  southern  tip  of  the  Kamchatka  peninsula  in  the  sea  of 
Okhotsk,  experiences  a  doubling  of  the  average  wind  stress 
and  a  179  W  m  ^  decrease  in  the  net  surface  heat  flux  due  to 
the  smoke  injection  (Table  2)  The  results  for  this  point  are 
show  n  in  Figures  4  and  5. 


Looking  first  at  the  atmospheric  forcing  (Figure  4).  we  see 
that  the  net  heat  flux  is  affected  by  the  smoke  as  early  as  July 
3,  and  that  an  exceptionally  strong  wind  event  lasting  several 
days  occurs  shortly  thereafter.  For  16  consecutive  days  during 
the  nuclear  winter  experiment,  from  July  1 1  to  26,  the  solar 
radiation  is  reduced  essentially  to  zero  by  the  thick  smoke 
cloud  originating  over  nearby  Siberia,  and  this  produces  a  net 
heat  loss  (Q„  <  0)  which  lasts  almost  to  the  end  of  the  month. 
By  the  end  of  the  month,  the  net  heat  flux  has  partially  re¬ 
turned  to  normal.  F.xcept  for  the  enormous  wind  storm  during 
the  first  Week,  the  wind  stress  is  quite  similar  in  the  two  exper¬ 
iments.  This  characteristic  of  a  strong  wind  event  early  in  the 
month,  followed  by  severe  reductions  in  solar  radiation  and 
subsequent  partial  recovery,  is  typical  of  the  atmospheric  forc¬ 
ing  during  nuclear  winter  at  other  locations  as  well. 

The  response  of  the  upper  ocean  thermal  structure  (E'igure 
5)  is  easily  understood  in  terms  of  the  above  forcing.  During 
the  control  experiment,  the  upper  ocean  responds  to  a  cycle  of 
typical  summer  wind  events  on  July  6,  13,  16,  21,  and  25.  The 
ocean’s  response  to  these  wind  events  is  seen  in  the  deepening 
and  coalescense  of  the  isotherms  representing  an  enhanced 
vertical  temperature  gradient  at  the  base  of  the  mixed  layer 
caused  by  the  wind-generated  downward  mixing  of  warm  sur¬ 
face  layer  water.  In  between  these  characteristically  gentle 
wind  events,  the  isotherms  tend  to  spread  out  in  the  vertical  as 
the  mixed  layer  shallows  and  warms  owing  to  the  net  surface 
heating.  The  evolution  during  the  nuclear  winter  experiment  is 
very  different.  The  wind  storm  on  July  5  deepens  and  cools  the 
mixed  layer,  and  the  subsequent  negative  (upward)  heat  flux 
and  relatively  normal  winds  produce  a  gradual  cooling  and 
slow  deepening  of  the  mixed  layer  during  the  remainder  of  the 
month.  In  spite  of  the  drastic  changes  in  the  atmospheric  forc¬ 
ing  during  the  nuclear  winter  experiment,  the  response  of  the 
ocean  is  confined  to  the  upper  60  m. 

Discussion 

In  this  study  we  have  examined  the  short-term  response  of 
the  upper  ocean  to  Ihe  simulated  forcing  of  a  hypothetical 
atmospheric  nuclear  winter.  The  atmospheric  forcing  was 
derived  from  fields  produced  by  a  two-level  atmospheric  gen¬ 
eral  circulation  model  in  both  a  control  and  a  nuclear  winter 
simulation.  The  ocean  was  represented  by  a  one-dimensional 
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mixed  layer  motiel  applied  al  six  difTerent  locations  in  the 
North  Atlantic  and  North  Pacitic  oceans.  Since  there  was  no 
negative  feedback  with  regard  to  the  heat  exchange  between 
the  ocean  and  atmosphere,  the  intensity  of  the  modeled  ocean 
response  may  be  regarded  as  an  upper  bound  (i.e..  the  actual 
response  may  be  weaker)  This  is  a  useful  thing  to  know. 


particularly  for  estimating  the  si/e  of  the  feedback  from  the 
ocean  and  its  potential  ctVcci  on  Ihc  atmosphere. 

Our  results  indicate  that  al  some  of  the  locations,  depend¬ 
ing  on  the  distance  upwind  from  the  point  in  question  to  the 
nearest  source  of  smoke,  the  upper  ocean  experiences  a  signih- 
cant  cooling  in  response  to  a  large  scale  nuclear  exchange. 


Ml  I  II  AI  M  I  I  M  .  Bull  I  RiI'IIKI  I4'' 


l  AHLl.  OitVcrencc  {Jeivvcfn  the  \K>del  Smiulaled  SST  al  the  Eind 

v'l  iht-  I'vso  fApcnmeniN  iC'ontrol  .Minus  Nuclear  Winter  = 
and  the  Dhsersed  Change  in  SST  During  a  Climatological  Seasonal 
C>cic  {Summer  .Minus  W  inter  ^ 

Point 

Location 

^Sfltv**ix 

i 

50  N.  150  \V 

:.xH 

4.0 

: 

50  N.  150  1 

5.0 

8.0 

's 

50  N.  .sO  W 

.V5 

40 

4 

.10  N.  150  \V 

1.0 

5.0 

5 

10  N,  150  W 

0.5 

20 

6 

.10  N.  125  F 

N.5 

12,0 

Ml  Icnipcnmirc  liill’ercnu"'  lire  in  degrees  Celsius,  with  elimiitologi- 
e.il  v.ilues  ffiim  Rohiiistui  [W76]  tor  the  North  Paeitic  Oeean  (points 
I.  -.4.  .s.  .iiiJ  M  and  from  KcynnUls  [1482]  for  the  North  Atlantic 
Ocean  ipoini  .'1. 


This  cooling  produces  a  decrease  in  the  SST  which  can  be  a 
si/eable  fraction  of  the  normal  seasonal  cycle  at  the  location 
111  question  (Table  .4|.  Thus  in  the  eastern  North  Pacific  Ocean 
along  150  W  the  decrease  in  SST  due  to  the  nuclear  exchange 
Is  approximately  .80  "i.  of  the  normal  seasonal  change.  At  the 
tnher  (hree  locations,  which  happen  to  lie  closer  and  down¬ 
wind  of  the  source  of  smoke  in  the  atmosphere  (Figure  1),  the 
decrease  in  SST  due  to  the  nuclear  exchange  is  about  75"'c.  of 
the  normal  seasonal  change.  With  the  exception  of  the  East 
China  Sea  location  discussed  earlier,  this  cooling  is  caused 
primarily  by  the  smoke  induced  decrease  in  solar  radiation. 
Additional  cooling  is  produced  by  enhanced  sensible  plus 
latent  heat  flux  and  entrainment,  but  the  combined  effect  of 
these  proces.ses  is  generally  smaller  than  the  effect  of  decreased 
solar  radiation  (Table  2i.  The  East  China  Sea  location  is 
unique  because  of  the  extremely  large  increase  in  the  sensible 
plus  latent  heat  flux  that  occurred  at  that  location  during  the 
nuclear  winter  experiment.  As  was  discussed  above,  this  very 
large  heat  loss  to  the  atmosphere  at  the  surface  is  considered 
to  be  somewhat  unrealistic  because  the  ocean  surface  temper¬ 
atures  were  kept  at  their  July  climatological  values  through¬ 
out  the  atmospheric  nuclear  winter  simulation.  A  more  realis¬ 
tic.  coupled  ocean-atmosphere  model  would  most  likely  pro¬ 
duce  smaller  air-sea  fluxes  and  thereby  result  in  a  weaker 


ocean  response  in  this  region.  Nevertheless,  the  relatively  large 
ocean  response  at  the  East  China  Sea  location  tentatively  sug¬ 
gests  that  locations  such  as  this,  having  rather  high  SSTs  ini¬ 
tially  and  lying  immediately  downwind  of  a  large  land  mass 
experiencing  a  nuclear  w  inter,  are  likely  to  experience  the  larg¬ 
est  SST  decreases  in  response  to  a  nuclear  exchange. 

Wc  have  presented  results  for  only  six  selected  geographical 
points,  an  admittedly  small  sample.  Whether  or  not  our  re¬ 
sults  arc  representative  of  other  oceanic  regions  can  be  esti¬ 
mated  to  some  extent  by  examining  the  geographical  distri¬ 
bution  of  the  smoke-induced  changes  in  the  mean  wind  stress 
and  surface  heating,  as  is  simulated  by  the  general  circulation 
model  (Figure  6).  Such  an  examination  indicates  that  both 
increased  wind  stress  and  decreased  surface  heating  arc  quite 
typical  of  the  response  to  the  smoke  for  the  northern  hemi¬ 
sphere  oceans;  regions  of  reduced  wind  stress  or  enhanced 
surface  heating  arc  generally  confined  to  land  surfaces  and  or 
the  southern  hemisphere.  Regions  of  particularly  intense  sur¬ 
face  cooling  are  found  in  the  Sea  of  Japan  and  in  the  Bay  of 
Bengal,  associated  with  enhanced  convective  fluxes  as  cold 
continental  air  drains  from  the  Tibetan  plateau.  Broad  bands 
of  enhanced  wind  stress  are  found  in  the  western  Pacific,  far 
eastern  Pacific,  and  western  Atlantic  oceans  at  10  N  latitude. 
They  reflect  a  westerly  acceleration  of  the  zonal  flow  by  the 
enhanced  Hadley  circulation  driven  by  a  smoke-induced  me¬ 
ridional  heating  gradient  [Coiev  et  ai.  1984]. 

At  the  present  time  i(  is  very  difficult  to  evaluate  the  signifi¬ 
cance  of  the  ocean's  short-term  response  to  the  kind  of  large- 
scale  nuclear  exchange  postulated  in  this  work,  and  it  is  vir¬ 
tually  impossible  to  estimate  from  the  present  results  what  the 
ocean's  longer-term  response  might  be  As  was  noted  earlier,  a 
large  fraction  of  the  SST  change  following  such  a  nuclear 
exchange  takes  place  during  a  very  short  period  of  time,  typi¬ 
cally  just  .several  days  (Figure  2).  This  is  quite  dilfcrent  from 
the  normal  seasonal  decrease  in  SST  from  summer  to  winter 
which  occurs  over  a  period  of  at  least  several  months.  Thus 
there  is  some  indication  that  1  month  after  the  smoke  is  intro¬ 
duced  into  the  atmosphere  the  SSTs  are  approximately  follow¬ 
ing  the  normal  annual  cycle  (i.e..  that  of  the  control  experi¬ 
ment)  only  at  a  lower  value  of  the  temperature.  However,  our 
experiments  are  too  short,  and  the  modeling  and  experimental 
uncertainties  are  too  great,  to  allow  us  to  infer  anything  about 
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changes  that  might  occur  in  the  ocean  on.  say,  the  seasonal 
time  scale.  Indeed,  some  of  our  results,  e  g.,  those  at  the  fast 
China  Sea  location,  suggest  that  it  may  be  necessary  in  the 
future  to  use  a  coupled  ocean-atmosphere  model  to  properly 
address  such  important  questions. 
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